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Figure 1. Fragment of the mural paint of the Xitle eruption and impact on Cuicuilco 
(Artist: Alejandro Linares-Garcia on display at Tlalpan Delegation, Mexico City). 


ABSTRACT. The Cuicuilco archaeological site in from the Xitle volcano. Dating the Xitle eruption and 
southern Basin of Mexico is covered by lava flows Cuicuilco abandonment has long been attempted. 
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Figure 2. Schematic map of Xitle lava flows, southern Basin of Mexico (Delgado-Granados et a/. 1998). Arrows indicate 
flow directions. Location of the Cuicuilco archaeological center is indicated (inset aerial view of Cuicuilco sector-A 
and the circular-cross section pyramid, surrounded by lava flows; Siebe 2000). 

The UNAM National University campus is built on the Xitle lava field. 


Contrasting results with radiocarbon dates around 
2000 and 1670 yr BP have been reported, with im- 
plications for the development of the Mesoamerican 
centers of Cuicuilco and Teotihuacan. Here, we 
analyze radiocarbon dates and paleomagnetic data 
for the Xitle lava flows. New age estimates for the 
eruption are determined from correlating full vector 
data with the geomagnetic secular variation refer- 
ence model. The revised archaeomagnetic data give 
ages correlating with the radiocarbon chronology, 
with a mean of 2086 cal yr BP and 95% confidence 
interval from 1995 to 2177 cal yr BP. Bootstrap analy- 
sis of the calibrated radiocarbon and archaeomag- 
netic dates gives mean dates and confidence inter- 
vals of 2041 and 1968-2041 cal yr BP and 2035 and 
1968-2073 cal yr BP, respectively. The interval es- 
timated of ~90 BC to ~AD 20 supports a possible link 
between the abandonment of Cuicuilco and the early 
development of Teotihuacan. 


KEYWORDS: Cuicuilco, Teotihuacan, Xitle volcano, 
Archaeomagnetism, Mesoamerica. 


RESUMEN. La zona arqueoldgica de Cuicuilco, en 
e/ sur de la cuenca de México, esta cubierta por flu- 
jos de lava del volcan Xitle. Se ha intentado la data- 
cidn de la erupcion y el abandono del centro de 
Cuicuilco aplicando diferentes métodos. Se han pro- 
puesto fechas contrastantes alrededor de 2000 y 
1670 anos AP, con implicaciones para el desarrollo 
de los centros urbanos mesoamericanos Cuicuilco 
y Teotihuacan. A continuacion, analizamos las fechas 
de radiocarbono y los datos paleomagnéticos para 
los flujos de Xitle. Se presentan nuevas estimacio- 
nes de la edad de la erupcidn usando datos del vec- 
tor completo con el modelo geomagnético de refe- 
rencia. Los datos paleomagnéticos revisados dan 
edades con una media de 2086 anos AP e intervalo 
de confianza del 95 % entre 1995-2177 anos AP. 
E/ analisis bootstrap de las edades radiocarbonicas 
y arqueomagnéticas proporciona edades medias e 
intervalos de confianza de 2041 y 1968-2041 anos 
AP y 2035 y 1968-2073 anos cal AP, respectiva- 
mente. El intervalo estimado de ~ 90 a.C. a 20 d.C. 
es compatible con una posible relacion entre el aban- 
dono de Cuicuilco y el desarrollo de Teotihuacan. 
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PALABRAS CLAVE: Cuicuilco, Teotihuacan, volcan 
Xitle, arqueomagnetismo, Mesoamerica. 


INTRODUCTION 


Xitle is a monogenetic cinder cone located in the 
southern Basin of Mexico, on the northern slope of 
the Ajusco volcanic complex. It forms part of the 
Chichinautzin monogenetic volcanic field developed 
in the arc front zone of the Trans-Mexican volcanic 
belt (Urrutia-Fucugauchi & Martin 1993). Xitle is one 
of the youngest volcanoes and its lava flows cover 
the Cuicuilco archaeological center (Figs. 1, 2). 
Cuicuilco represents an early urban settlement in the 
Basin of Mexico (Heizer & Bennyhoff 1958). The 
eruption of Xitle volcano has been linked to the aban- 
donment of the site, which was covered by tephra 
and lava flows. Beginning with the initial excavation 
projects, determining the age of the eruption has 
attracted considerable interest. Arnold and Libby 
(1951) reported a date of 2422 + 250 yr BP, as part 
of the first date set obtained with the radiocarbon 
method (Libby 1955). Since then, more than 30 ra- 
diocarbon dates for samples from the Xitle lavas and 
Cuicuilco archaeological site have been reported, 
which present a wide multimodal ~4000 to ~1500 yr 
BP distribution (Urrutia-Fucugauchi 1996). 

The first date of Arnold and Libby (1951) was ob- 
tained from a soil sample recovered beneath the 
lava. Subsequent studies have dated samples from 
soils and charcol from outcrops and from archaeo- 
logical excavations in different sectors of the lava 
field and in the Cuicuilco site. Dates around 2000-— 
1960 yr BP and younger dates around 1530-1630 
yr BP have been proposed for the Xitle eruption (e.g. 
Fergusson & Libby 1963; Cordova et al. 1994; Urru- 
tia-Fucugauchi 1996; Siebe 2000; Gonzalez et al. 
2000). Older or younger dates present archaeologi- 
cal implications for the abandonment of the site, 
population migration and relationships to the devel- 
opment of the Teotihuacan urban center (Nichols 
2016). Discussion on the age of the eruption has fo- 
cused on the stratigraphic context and nature of the 
samples dated, with arguments for and against its 
temporal relation to the tephras and lavas (Cordova 
et al. 1994; Gonzalez et a/. 2000). In this study we 
present an analysis using paleomagnetic data re- 
trieved from the lava flows. For the archaeomagnetic 
dating we use full vector data with the remanent mag- 
netization directions and paleointensities and the re- 
cently developed archaeomagnetic reference curve. 
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ARCHAEOMAGNETIC AND 
RADIOCARBON DATING 


The basaltic flows of the Xitle eruption cover an ex- 
tensive > 80 km? area on the Ajusco volcano slope 
and into the adjacent plain, where the Cuicuilco 
center was built and nowadays is occupied by the 
urban development of Mexico City (Fig. 2). The flows 
have long been studied, with several distinct flow and 
tephra units mapped (Delgado-Granados et al. 
1998). Paleomagnetic studies have been conducted, 
reporting remanent magnetization directions and 
paleointensities from different sites in the lava field 
(e.g. Nagata et a/. 1965; Urrutia-Fucugauchi 1996; 
Morales et al. 2001, 2006; Alva-Valdivia 2005). 

The paleomagnetic directions present northward 
declinations and downward inclinations with small 
within-site angular dispersion. The mean-site incli- 
nations mainly cluster around the dipolar direction, 
with some showing shallower inclinations (Urrutia- 
Fucugauchi 1996). The overall mean direction de- 
termined for 26 sites is Dec = 359.8°, Inc = 32.8°, k 
= 167 anda,, = 2.2°. The mean direction for the shal- 
low inclination group is N = 7, Dec = 1.5°, Inc = 26.7°, 
k = 338 and a,, = 3.3°. The mean direction for the 
dipolar inclination group is N = 19, Dec = 359.02, Inc 
= 35.1°, k= 247 anda,, = 2.1°. Paleointensities were 
determined for some sites, using the Thellier dou- 
ble-heating and Shaw methods, giving a mean of N 
= 20, 60.5 + 9.2 uT (a subset group gives a mean N 
= 13, 56.5 + 6.2 uT). Alva-Valdivia (2005) reported 
data for 10 sites, with a mean direction of Dec = 1.1°, 
Inc = 34.1°, k = 531 and a,, = 2.1°. In this study, a 
larger set of paleointensities determined with the 
Thellier method was obtained. The paleointensities 
range from 48.6 to 73.9uT, with mean of N = 9, 59.9 
+ 7.7uT. For this study, the site mean directions and 
paleointensities reported in Urrutia-Fucugauchi 
(1996) and Alva-Valdivia (2005) are referred to the 
archaeomagnetic reference curve by Pavon-Car- 
rasco et al. (2014). Inclusion of full vector data im- 
proves analysis, considering three variables with 
higher-resolution age estimates. 

The radiocarbon dates range from 4690 + 70 to 
1430 + 200 yr BP, with four groupings around 4000, 
2400, 2000 and 1500 yr BP. Production of radioac- 
tive carbon varies with the strength of the geomag- 
netic field, which changes through time (Ramsey et 
al. 2006). Correcting the dates using the dendro- 
chronologic reference curve provides calendar years 
ranging from about 3500 BC to AD 640, which spans 
a rather long interval. Radiocarbon dates ranging 
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Figure 3. Archaeomagnetic dating of Xitle lavas. Diagrams show in red the master reference archaeomagnetic curves with un- 
certainty ranges for declination, inclination, and intensity (Pavén-Carrasco et a/. 2014 geomagnetic model); and in blue the Xitle 
paleomagnetic data. Green lines give the cutoff values at a given probability level. Probability density functions from the secular 
variation curve analysis are shown in the lower diagrams. The combined analysis of full vector paleomagnetic data is given in 
the summary diagram with the reference geographic map. Results for sites XT-7 (Urrutia-Fucugauchi 1996, see above) and 6 
(Alva-Valdivia 2005, see below). 


from 1950 to 2100 yr BP result in calendar years were not involved in the center expansion (Nichols 
around 100 BC to AD 60, which correlate with the 2016; Cordova et al. 1994; Gonzalez et a/. 2000). 
expansion of Teotihuacan. Younger radiocarbon _ For this study, radiocarbon dates are recalibrated to 
ages of about 1670 yr BP (Siebe 2000) imply that the dendrochronologic curve and re-analyzed using 
Cuicuilco inhabitants, if migrating to Teotihuacan, bootstrap statistics. 
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Figure 4. Distribution of radiometric (a) and archaeomagnetic dates (b). Mean dates and confidence 
limits derived from bootstrap analysis. Insets show Gaussian function fits (see text). 
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ANALYSIS AND DISCUSSION 


Archaeomagnetic dating relies on the spatio-tempo- 
ral variations of the geomagnetic field and the record 
of remanent magnetization in volcanic and archaeo- 
logical materials. Remanent magnetization, direc- 
tional and paleointensity data, which record the geo- 
magnetic field at time of magnetization acquisition, 
are correlated to a reference curve independently 
dated from the directional and intensity changes of 
the geomagnetic field. Resolution depends on the 
fidelity of the remanence recording and precision of 
the reference geomagnetic secular variation curve. 
For this study, the paleomagnetic data for the Xitle 
lavas have small angular dispersion at the within- and 
between-site levels. The mean directions determined 
in Urrutia-Fucugauchi (1996) and Alva-Valdivia 
(2005) are statistically indistinguishable, with small 
2.1° a,, overlapping confidence cones. Site-mean di- 
rections show northward declinations and inclina- 
tions close to the dipolar inclination, which are re- 
ferred to the geomagnetic reference SHA.DIF.14k 
model of Pavén-Carrasco et al. (2014). We consider 
full vector data, adding the paleointensities in the cor- 
relation. The paleointensities present higher disper- 
sion at the between site level, which has been ana- 
lyzed, incorporating determinations with cooling rate 
corrections (Morales et al. 2006; Alva-Valdivia 2005). 

For this study, the archaeomagnetic correlations 
are performed at site level (Fig. 3), to take into ac- 
count the dispersion. The results show internal con- 
sistency, with small standard deviations. Assuming 
a Gaussian probability distribution, mean is 2068.5 
cal yr BP with sigma of 143 and 95% confidence in- 
terval from 1995.4 to 2177.4 cal yr BP (Fig. 4). Analy- 
sis of radiocarbon dates in the 1200-2600 cal yr in- 
terval gives 12 dates with mean of 1916.8 cal yr BP 
and sigma of 227, with 95% confidence interval from 
1772.2 to 2061.3 cal yr BP (Fig. 4). The mean ra- 
diocarbon date is younger than the archaeomagnetic 
date. 

To analyze the difference in age estimates, we 
carried out a bootstrap analysis to constrain the best 
date estimates. Confidence limits are determined for 
each date and used to generate random series of 
100 vectors. For the radiometric dates, the resulting 
matrix then has 12 x 100 elements and for the ar- 
chaeomagnetic dates it has 9 x 100 elements. Vec- 
tors are fitted using splines and mean and confi- 
dence intervals calculated. The resulting mean for 
the archaeomagnetic dates is 2035 cal yr BP, witha 
95% confidence interval from 1968 to 2073 cal yr BP. 
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The corresponding mean for the radiometric dates 
is 2041 cal yr BP, with 95% confidence interval from 
1968 to 2041 cal yr BP (Fig. 4). 

These estimates can be refined by deleting data 
points that fall off from the confidence interval in the 
probability density distributions. In this analysis, one 
archaeomagnetic and two radiometric dates are 
deleted. The resulting mean for the archaeomagnetic 
dates is 2101 cal yr BP with confidence interval of 
2065-2137 cal yr BP. For the radiocarbon dates the 
mean is 2061 cal yr BP with a confidence interval of 
1985-2132 cal yr BP. 


CONCLUSIONS 


The age of the Xitle eruption is determined from 
correlating paleomagnetic full vector data with the 
recently constructed geomagnetic secular variation 
reference model. The revised archaeomagnetic 
dates have a mean age of 2086 cal yr BP with 95% 
confidence interval from 1995 to 2177 cal yr BP. 
Bootstrap analysis of the calibrated radiocarbon and 
archaeomagnetic dates gives mean ages and con- 
fidence intervals of 2041 cal yr BP and 1968-2041 
cal yr BP and 2035 cal yr BP and 1968-2073 cal yr 
BP, respectively. These estimates are internally con- 
sistent, with overlapping calendar intervals of 91 BC 
to AD 18 and 83 BC to 18 AD, respectively. Results 
support a possible link between the Xitle eruption and 
collapse of Cuicuilco center and the early develop- 
ment of Teotinuacan. 
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